The present contribution reports some preliminary results obtained applying a simple finite element formulation, developed for discretizing the partial differential equations of motion of a novel beam model. The theoretical model we are dealing with is geometrically exact, with some peculiarities in comparison with other existing models. In order to study its behavior, some numerical investigations have already been performed through finite difference schemes and other methods and are reported in previous contributions. Those computations have enlightened that the model under analysis turns out to be quite hard to handle numerically, especially in dynamics. Hence, we developed ad hoc the total-lagrangian finite-element formulation we report here. The main differences between the theoretical model and its numerical formulation rely on the fact that in the latter the absolute value of the shear angle is assumed to remain much smaller than unity, and strains are piecewise constant along the beam. The first assumption, which actually simplifies equations, has been taken on the basis of results from previous integrations, mainly through finite difference schemes, which clearly showed that, while other strains can achieve large values in their range of admissibility, shear angle actually remains small. The second assumption led us to define a two-nodes constant-strain finite element, with a fast convergence, in terms of number of elements versus solution accuracy. Although, at the present stage of this ongoing research, we have only early results from finite elements, they appear encouraging and start to shed new light on the behavior of the beam model under analysis.
Explicit finite difference schemes for solving the one – dimensional tsunami wave propagation equations
